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ABSTRACT 

The Meshless Local Petrov-Galerkin (MLPG) 
method is one of the recently developed element-free 
methods. The method is convenient and can produce 
accurate results with continuous secondary variables, 
but is more computationally expensive than the finite 
element method. To overcome this disadvantage, a 
simple Heaviside test function is chosen. The 
computational effort is significantly reduced by 
eliminating the domain integral for the axisymmetric 
potential problems and by simplifying the domain 
integral for the axisymmetric elasticity problems. The 
method is evaluated through several patch tests for 
axisymmetric problems and example problems for 
which the exact solutions are available. The present 
method yielded very accurate solutions. The 
sensitivity of several parameters of the method is also 
studied. 

INTRODUCTION 

The Meshless Local Petrov-Galerkin (MLPG) 
method is a promising numerical method to analyze 
potential and elasticity problems and is shown to yield 
accurate results [1-6]. In this method, a set of 
arbitrarily distributed nodes is used to interpolate the 
field variables. The method does not use either 
‘elements’ or a back-ground mesh for integration and 
hence is truly meshless. The method, however, 
appears to be more computationally expensive than 
the finite element method. One of the reasons for the 
high computational cost is that accurate numerical 
integration is required to integrate the weak form of 
the governing equations. Thus, if the domain 
integrations are eliminated or simplified, the method 
can be made efficient. A choice of a “Heaviside” test 
function [5, 6] leads to the elimination of the domain 
integral in the weak form of the axisymmetric 
potential problems and substantially simplifies the 
integral for axisymmetric elasticity problems. The 
purpose of this paper is to present such a method that 


utilizes the Heaviside test function and evaluate its 
effectiveness for axisymmetric problems. 

The outline of the paper is as follows. First, a brief 
overview of MLPG method is presented for 
axisymmetric potential and elasticity problems. The 
MLPG formulation with Heaviside function is 
presented, and various issues related to the use of a 
Heaviside test function are studied. Finally the 
effectiveness of the method is demonstrated using 
several numerical problems. 

OV ERMEW OF MLPG METHOD 

In this section, several basic concepts that are used 
in the MLPG method are briefly reviewed. The 
development of the weak forms of the solution for 
potential and elasticity problems is first presented. 
Next, the choice of the trial function for the primaiy 
variables using the moving least squares method is 
reviewed. In the classical element-free Galerkin 
methods, the test functions are chosen from the same 
space as the trial functions. In the current MLPG 
method, the test functions are chosen to be different 
from the trial functions, and the choice of the test 
functions is discussed next. 

Potential Problem 

Consider Poisson’s equation for an axisymmetric 
problem bounded by a toroidal domain with its cross 
section defined by Q as shown in Figure 1, 

V 2 u = g in Q (1) 

with boundary conditions 

u = u on T u and q-q on F q (2) 

where T =T U + r q and q = du / dn . The Laplacian in 
the cylindrical coordinate system is 
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and for the axisymmetric problem, Eq. (3) reduces to 


d 2 1 d d 

T ^ 1 7 

dr" r dr dz~ 


The solution for Eq. (1) is sought in a weighted 
residual manner as 


where F s is the boundary of the compact support 
domain and T is O (see Fig. 1) and q is 
the prescribed flux. 

Elasticity Problem 

Consider a problem of axisymmetric elasticity in a 
toroidal body with its cross section defined by domain 
Q (Fig. 1), bounded by T. The equations of 
equilibrium are 


2n J (V 2 t/ -g)-v-r-dr-dz = 0 (5) 

Q 

where u is the approximate solution of the boundary- 
value problems and v is a weight function. 

Using the divergence theorem, one can rewrite Eq. 
(5) (dropping the constant 2 n) as 


-I 
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dr dr dz dz 


■ r ■ dr ■ dz 
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dr 

dr 
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( 8 ) 


in Q where <7 r ,<J g ,(J_ and l r2 are the stresses in 
cylindrical coordinates ( r,9,z ) and f. ,f z are the body 
forces in radial (r) and vertical (z) directions, 

respectively. The boundary conditions on F can be 
written as 


Uj = M,. on r u 


(9) 


where n, and n, are the direction cosines of the normal 
to the domain boundary along the r- and z-directions, 
respectively. Eq. (6) is the weak form of the governing 
equation for potential problems. 

In the literature, the interpolations for the primary 
variables u are customarily termed as trial functions 
and the choice of the weight functions v as the test 
functions. In the MLPG method, these choices are 
made from different spaces. The test functions are 
chosen such that they are non-zero over an arbitrary 
shape of domain Q s that is known as compact support 
domain. In general, the domain is chosen to be a well- 
defined geometric shape such as a circle, a square, a 
rectangle or an ellipse. (In this paper, a compact 
support in the shape of circle is used.) As such, the 
integrals over Q in Eq. (6) reduce to integrals over Q s . 
Eq. (6) can be written as 


I 

Q s ' 

I 


du dv du dv 
dr dr dz dz 


r ■ dr ■ dz + 


du du | „ 

— n - 1 n_ | r- v-di + 

dr dz 


(7) 


J q ■ r- v ■ dT - J* g ■ v ■ r ■ dr -dz = 0 


&ij-n j = t i onT q (10) 

where rij are the direction cosines to the normal of the 
boundary, and U and t are prescribed displacements 
and tractions, respectively. The weighted residual 
form for axisymmetric elasticity problems can be 
obtained using 

2 A ^ + -(a r - a e ) + f.) ■ v r ■ r ■ dr ■ dz = 0 
J Q or oz r 

(ii) 


In [ (Aa + + ~z r _ + /_) • v. • r ■ dr -dz -0 (12) 

•j dr dz r 

where v r and v, are weight functions for r-and z- 
directions, respectively. Using the divergence 
theorem, Eqs. (11) and (12) can be rewritten to obtain 
the weak form of the Eq.(8) as 


f / ov r 3v cr„ . , , 

• I Mr ■ -r - + ^ + — ■ v.) • r ■ dr ■ dz - 

•'L2 rlr Qz y 


dv r 
dr 
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(13) 


2 

American Institute of Aeronautics and Astronautics 



r ov Ov 
- [ (r : -+(7 — -)-r-dr-dz + 

rz dr z dz 

[ f -v -r-dr-dz+ (1^) 

Jsr 2 - 

J r (*W», +VM : )-v z -r-dT =° 

By choosing the test functions v, and v- to be non- 
zero only inside a compact support domain Q s , the 
weak form (Eqs. (13) and (14)) can be rewritten as 


j n ( ',<r r ■ ^ ^ + — ■ F ) f dr • dz + 

Or dz r 

L (<FA + Fz«-) ' F ' >' ' «T + (15) 

Jl S 
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where t are prescribed tractions on F . 

Trial functions using Moving Least Squares 
Interpolation 

The moving least squares (MLS) approximations 
have reasonably high accuracy and can be extended to 
/7-dimensional problems, and hence are widely used to 
construct the trial functions in meshless methods [1- 
6]. This interpolation is used to represent the trial 
functions of the primary variables at a set of randomly 
located nodes in the domain. The interpolation 
function i/\x) of MLS method is defined as 


u h (x) = p r (x)-a(x) VxeQ (17) 

where p T (x) is the monomial basis function of order 
m. For axisymmetric problems, monomial basis 
functions are chosen as 


P T (x) 


[ 1 ,r,z\ 

[ 1 ,r,z,r 2 ,rz,z 2 ] 

[ 1 ,r,z,r 2 ,rz,z 2 ,r 3 ,r 2 z,rz 2 ,z 3 ] 


(m = 3) 
(m = 6) 
(m = 1 0) 
(18) 


where m = 3, 6, and 10 are complete linear, quadratic 
and cubic basis functions, respectively. 


The coefficient vector a(x) in Eq. (17) can be 
obtained by minimizing the function J(\) that is 
evaluated at n discrete nodes 


j (x) = Tj W, (x) • [p r (x,. ) • a(x) - uf , (19) 

!= 1 


where U is a vector of the fictitious nodal value of u, 
and vv,(x) is a weight function. (Figure 2 illustrates the 
MLS scheme for a 1-D example.) In this paper, two 
spline functions, one with C 1 and another with C 2 
continuity, are used as weight functions. The weight 
functions chosen are: 


C 1 : 

w*(x) 


- 3 A- 2 + 2 A- 3 0<p k <\ 

0 A > 1 


(20) 


C 2 : 

AW 


-10a 3 + 15A 4 - 6A 5 0 < a < 1 (21) 
0 A > 1 


where p k = (d k / l k ) is the normalized distance, 
d k = ||x - x k | is the distance between center of the 

sub-domain of influence \/ ( (r k , z k ) and a general point 
x, and l k is the radius of the domain of influence 
(extent of the trial function) at node k (see Figure 1). 

By minimizing Eq.(19), then substituting the result 
into Eq.(17), the MLS interpolation function (see Fig. 
2 for a one-dimensional example) is obtained as [1-6] 


u\x) = Y J (pi^)-u i ( 22 ) 

i=i 

where n is the number of nodes in the domain, u j is 
the fictitious nodal value of u at node i, and ^r(x) is 
the shape function of node i, and is defined as 

m 

</,(*) = (23) 

r = i 

with 


A(x) = X w /( x ) ' P r (x,) -p(x,) (24) 

1=1 

B(x) = [tVj (x) • p(Xj), w 2 (x) ■ p(x 2 ), ...w n (x) • p(x n )] 

(25) 
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Note that in the MLS approximation u. evaluated 
using Eq. (22) is not equal to the fictitious nodal vaule 
U j (also see Figure 2). The partial derivatives of 

(Pj{x) are needed in the weak form (Eqs.(7), (15) and 
(16)). Since both matrix A(x) and B(x) are functions 
of x, the partial derivatives of ^>,(x) have 
complicated form as [2] 

m 

<P iJc = Il/f, x (A B), . /,,(A B. iA,B).| (26) 

7=1 

where 

A-^CA^^-AV' (27) 

Note that (p f , in Eq. (26) represent (p i r and (p. f z . 


- f — • v„ • r ■ dr ■ dz + f /' • v„ • r-dr-dz + 

Jn s r Jo, r r 

j r (an r + T r n z )-Vr-r-dr =0 

[ f-v ■ r-dr-dz + 

J n 5 z - 

j, ( T ,~ n , + °VO • v z -r-dY = ° 

In this paper, the MLPG algorithm using the 
Heaviside test function is applied to various 
axisymmetric problems to study salient features of the 
method. The sensitivity of the MLPG solution to 
various parameters such as the size of the compact 
support domains, the basis function, etc., will be 
presented to demonstrate the algorithm. 

EXAMPLES 


Test function using Heaviside Function 

The next step in the MLPG method is to choose the 
test function. In the current implementation, a 
Heaviside function test function is chosen as the test 
function 


v(x) 


Jc x e Q s 

[o xgn 


(28) 


where c is an arbitrary constant (c = 1 is used in this 
study). Using this choice, the partial derivatives of the 
test function dv / dr and dv / dz are identically zero, 
hence the domain integral involved in Eqs. (7) and 
(16) are identically zero, and Eq. (15) for elasticity 
problems is also considerably simplified. The 
elimination of domain integrals over f2 s in Equations 
(7) and (16) considerably reduces the computing 
effort and hence makes this method computationally 
efficient. 

For potential problems, Eq. (7) is simplified as 


\ 


du du 

— n,. h n. 

dr r dz - 


■ r- v-dY+ J q ■ r- v-dY - 


j* g-v-r-dr-dz = 0 


(29) 


and for elasticity problems, Equations (15) and (16) 
reduce to 


To evaluate the current MLPG algorithm for 
potential and elasticity problems, several patch test 
problems with exact solutions are considered. The 
sensitivity of the MLPG solution to the basis functions 
and the two weight functions w k (R ) are studied. Since 
all of the potential and elasticity problems in this 
study are linear problems, the dimensions of the patch 
test models are multiples of a constant a. Various 
normalized values for the radii of domain of the 
influence, ( Ifa ), and the radii of the support domains, 
(/„/«), are used to study their influence on the accuracy 
of the results. 

The results of the MLPG method are compared to 
exact solutions. The error norm (||e*/H 2 ) is used to 
evaluate the effects of various parameters, and is 
given by 


IML = j^{( U j- U exact) /u exact } 2 j M ’ ( 32 ) 

where M is the total number of randomly distributed 
internal points in the domain at which the numerical 
solution is evaluated and compared to the exact 
solution, u exact . Note that these internal points are 
independent points and are not associated with the 
nodes used in the MLPG models. A value of M = 50 
is used for all models. 

Potential Problem 
Patch tests 

Two patch tests involving Laplace- ( V 2 m = 0 ), and 
Poisson- ( V 2 m = constant ) equation problems are 
considered. Figure 3 shows two models used in the 
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patch tests of the Laplace and Poisson equation 
problems. These examples are analyzed using the 
MLPG algorithm for various combinations of 
prescribed boundary conditions involving u and q. The 
Laplace equation over a triangular domain with two 
mixed boundary conditions is studied using a 
quadratic basis function. An exact solution for this 
problem is u = r 1 - 2z 2 + 3z . The radius of the 
domain of influence l k is set to 4 a. Two different 
choices for l a are considered. In the first choice, the 
distance from each internal node to the closest 
boundary is computed. These distances are used as the 
values of I a for that node. Using this definition, the 
support domain radius may be different for each node. 
The MLPG algorithm with different values of l a for 
each node recovered the exact solution for the Laplace 
equation problems. In the second choice, the smallest 
value of all the /„ values from the first choice is used 
for all the nodes. Once again, the exact solution was 
recovered to machine accuracy. In these studies, the 
effect of the radius of the domain of influence, l k , was 
also evaluated by varying the radius from 4 a to 10a. 
The accuracy of the results is insensitive to the radius 
of the domain of influence, I k . 

The sensitivity of the MLPG solution is further 
studied by varying the support domain radius and 
the error norm is shown in Figure 4. All of the nodes 
in the model were set to have the same value of I a . 
The results show that, for the Laplace equation, the 
two weight functions considered with C 1 and C 2 
continuity (Eqs.(20) and (21)) gave excellent solutions 
for both boundary value problems and for a variety of 
the normalized support domain radius, IJa. 


where /q and r 2 are the inner and outer radii, 
respectively. The heat flux, q, is defined as 


<7 


( t x -t 2 ) 

r - ln(/- 2 //-) 


(34) 


The MLPG model with 25 randomly distributed 
nodes and the boundary conditions are also shown in 
Figure 6. The exact solution for the temperature is a 
natural logarithmic function, and since the natural 
logarithmic function cannot be exactly represented by 
a polynomial, a study of the use of various basis 
functions (polynomial functions) is undertaken. Three 
polynomial basis functions, linear, quadratic, and 
cubic, are considered (Eq. (18)). Note that the higher 
order basis function requires substantially more 
computational effort than the lower order basis 
function. Figure 7 shows the error norm for various 
basis functions with a C 1 weight function. The support 
domain radius, l 0 , is set equal to 0.35a for the nodes 
on the boundary and equal to the distance to the 
closest boundary for the internal nodes. Accurate 
solutions are obtained for all cases studied. Figure 7 
shows that the higher order basis functions yielded 
more accurate results than the lower order basis 
functions. 


Elasticity Problems 


For all the elasticity problems studied, a Young’s 
modulus, 77 = 100 GPa, and a Poisson’s ratio, V =0.3 
are used. 


Patch Test 


The next problem considered involves a Poisson’s 
equation, V 2 m = 8 , over the L-shaped domain shown 
in Fig. 3. An exact solution for this problem is 
u = r 2 + 2z 2 + 3z - 2 . The convergence of the 
solution for various values of l 0 and for different 
boundary conditions is presented in Figure 5. Once 
again, the exact solution is recovered by MLPG 
algorithm for all idealizations demonstrating that the 
present MLPG algorithm is very accurate. 

Heat Transfer Example 


The axisymmetric elastic MLPG method was 
evaluated for patch test problems with the 
displacements ( u r =c x -r ; u_ = c, • z ) of a circular 

cylinder with different combinations of boundary 
conditions. Since exact solutions for these problems 
can be represented by polynomial functions, the 
MLPG method recovered the solutions to machine 
accuracy for all weight and basis functions considered. 
As in the cases of the potential problems, the accuracy 
of the results is insensitive to the radius of the domain 
of influence, l k . 


Figure 6 shows a heat transfer problem involving 
radial heat flow in a hollow circular cylinder with 
prescribed constant temperatures 7) and 77 at the inner 
and outer surfaces, respectively. The exact solution for 
the temperature at radius r is defined as 


T = T { - (Tj - J 2 ) • 


ln(r / r x ) 

InfeAi) 


(33) 


Lame’s Cylinder 

A hollow cylinder subjected to constant internal (p{) 
and external (p 2 ) pressures (Lame’s Problem) is 
considered (see Figure 8). For the problem analyzed, 
the inner radius r t and the outer radius r 2 are 3 a and 
6 a, respectively, where a is an arbitrary constant. The 
nodal arrangement is shown in Figure 8b. The exact 
solution of Lame’s problem is 
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cr. = Cj - CJr 1 ;a e = C l + C 2 /r 2 ;a z =2-v-C l 


1 + v 


C, 


= - 


Cj • (1 - 2v) • r + — 


m_=0; (35) 


where 


uniformly rotating long cylinder. Both inner and outer 
surfaces are stress free, while the disk is subjected to 
an inertial force per unit volume that equals p ■ o ) 2 ■ r , 

where p is the density of the disk material and a> is the 
angular velocity. The exact solution of the problem is 
given by (see pages 335-337, Ref. [7]) 


_ Pr>\ 2 -p 2 -r 2 . r _ P\ ~ Pi 2 2 

'-'l __2 _.2 ’ ^2 2 2 M '2 


2 2 
r 2 -r x 


r-r 

'2 'l 


(3-2 -k) 2 2 2 

cr,. = y9&r(/| + r 2 ■ 

r 8 • (1 - v) 


r 2 ) 


Figure 9 shows the variation of the error norm for 
three cases using a linear basis function. In the first 
two cases, the trial functions are formed using C 1 and 
C 2 continuity weight functions. The Heaviside test 
function is used in both cases. In the third case, the 
linear basis function and C 1 continuity weight function 
are used to form the trial function, and a C 1 spline test 
function is chosen as in reference [4]. All of the nodes 
in the model were chosen to have the same support 
domain radius, l 0 . As seen in the figure, the error norm 
varies substantially as the support domain size 
changes for all of the cases. For the third case that 
utilizing the C 1 spline test function, the error norm 
decreases monotonically as the normalized support 
domain radius increases. However, for the problems 
that used Heaviside test function, the error norm does 
not decrease monotonically. The error norm for the 
cases using the Heaviside test function reaches a 
minimum value near IJa = 0.6, then begins increasing 
again as the support domain radius increases further. 
The accuracy decreases for values of ( IJa ) > 0.6. 
This type of behavior was noted only for the results 
obtained using the Heaviside test function. Figure 9 
also shows that the error norm of the solutions 
obtained using the C'-weight function is lower than 
that with the C 2 -weight function for linear basis 
function. 

The effects of the basis function are also evaluated. 
The error norm obtained with each of the basis 
functions is presented in Figure 10. The exact 
solution for the displacement in Lame’s problem 
cannot be represented by any of the three - linear, 
quadratic and cubic - polynomial basis functions. 
Unlike for most problems studied, the linear basis 
function yielded the most accurate results, while the 
quadratic and cubic basis functions produced similar 
results. 

Rotating Disk 

The next problem studied is that of a hollow disk 
rotating about the z-axis with a radial velocity of m 
rad/sec as shown in Figure 11. In this problem, both 
top and bottom surfaces are constrained in the z- 
direction (“plane strain” condition) to simulate a 


(7a 


(3-2-vQ 

8-(l-v) 


• pa> 2 (r 2 + r 2 2 + 



a, = v- (a r + a 0 ) ; «_ = 0 


1 + 2 ' V r 2 ) 
(3-2 -v) 


u, = r ■ s Q 


(l + v)-(3-2v) 2 

pco 

8 E 


, 2 2 , r 2 -r 2 (l + 2v) 3 

r ■ (r x + r 2 ) + J - - r 


r (3 - 2v) 


v ■ (1 + v) ■ (3 - 2v) 2 

pa>~ 

8(1 -v)E 


P , , 2 ^ r A ■>; 


r-( r t +r 2 )- 


(36) 


where E is Young’s modulus of the material and v is 
Poisson’s ratio. 


The three models in Figure 11 are used to study the 
convergence of the solution to nodal refinement. In 
Model- 1, seven nodal points are placed in the r 
direction. In Model-2, the distance between adjacent 
nodes is half of the distance between adjacent nodes in 
Model- 1 (A/2). In Model-3, the distance between 
adjacent nodes is a quarter of the distance in Model- 1 
(A/4). The number of nodes in the z direction is kept 
the same (equal to 5) for all of the models. 

Figure 12 shows the effect of the nodal arrangement 
on the error norm. Model-3 has the most nodal points 
and produced, as expected, the most accurate solution. 
However, the analysis using Model-3 requires 
substantially more CPU time than that for Model- 1. 
The effect of the basis functions is shown in Figure 13 
using Model-2 with the C 1 weight function. The cubic 
basis function yielded a more accurate result than 
either the linear or the quadratic basis function. 
However, as shown in Figure 13, the analysis using 
the linear basis function is only marginally less 
accurate than the cubic basis function and the linear 
basis requires considerably less computational effort 
than the analysis using the higher order basis function. 
Thus, the use of the linear basis function is preferred 
for this problem. 

Figure 14 shows the effect of the support domain 
size on the error norm using Model-2 and a cubic 
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basis function. The results, once again, show that the 
error norm pass through a minimum value and then 
begins to rise as the support domain radius increases 
past a value l fa = 0.18. 

CONCLUDING REMARKS 

A meshless Local Petrov-Galerkin (MLPG) method 
that uses a Heaviside test function is presented and 
used to analyze potential and elasticity problems in 
axisymmetric domains. By using a Heaviside test 
function for the axisymmetric potential problems, the 
domain integral in the weak form is eliminated and a 
line integral is used instead. This substantially reduces 
the computation effort to construct the ‘stiffness’ 
matrix and hence is computationally efficient 
compared to the conventional MLPG method. In the 
axisymmetric elasticity problems, since the domain 
integral is simplified but is not eliminated, the 
reduction of the computational effort is limited. The 
MLPG algorithm using a Heaviside test function is 
implemented and its efficiency and accuracy are 
studied with respect to three basis functions, two 
different weight functions, and different sizes of the 
local support domain. 

Potential patch test problems involving Laplace and 
Poisson equations are used to evaluate the efficiencies 
of the algorithm. To evaluate the effectiveness of the 
MLPG algorithm for elasticity problems, several 
simple patch test problems, a Lame cylinder problem 
and a rotating disk problem, are studied. 

The C 1 weight function yielded more accurate 
results than the C 2 functions for the patch test 
problems studied. Since the analysis using the C 1 
function requires less computing time than that using 
the C 2 functions, the C 1 function is recommended. 
Unlike the MLPG method that uses a local spline test 
function (Ref. [4]), the error norm for the MLPG 
method with the Heaviside test function passes 
through a minimum value as the support domain 
radius, l a , increases. For the problems studied, the 
variation of the radius of the domain of influence, 4, 
does not show a significant effect on the accuracy of 
the results. For the problems studied, the use of the 
cubic basis function did not always yield the most 
accurate results and substantially increases the 
computing time compared to the use of the linear basis 
function. In general, increasing the number of the 
nodes in the model increases the accuracy of the 
result. 

In summary, the MLPG algorithm with a Heaviside 
test function yielded accurate solutions for the 
potential and elasticity patch tests as well as the other 
problems studied. A simple C 1 weight function, and a 


simple linear basis function are recommended to 

obtain accurate results. 
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Figure 1 Typical nodal pattern for an axisymmetric potential problem 



Figure 2 Approximation for u using MLS scheme for a 1-D example 
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(b) L-shaped domain 


Figure 3 Patch tests models for Laplace and 
Poisson’s equations 



(a) Boundary value problem -1 



Normalized support domain, /</a 
(b) Boundary value problem -2 


Figure 4 Laplace problem with a triangular domain - Convergence of the solution for various 
values of l 0 ( V 2 n = 0 ; Quadratic basis; 4= 5a) 
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Normalized support domain, IJa 
(a) Mixed boundary value problem -1 



(b) Mixed boundary value problem -2 


Figure 5 Poisson’s problem with an L-shaped domain - Convergence of the solution for 
various values of l 0 ( V 2 u = 8 ; Quadratic basis; 4= 4 a) 



Figure 6 Heat transfer in a hollow circular cylinder 
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Figure 7 Heat transfer in a hollow circular cylinder - 
Effect of basis functions (C 1 weight function) 


(a) Lame’s cylinder 



fb) Meshless model 
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Figure 8 Hollow cylinder subjected to internal and 
external pressures (Lame’s Problem) 
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Figure 9 Lame’s cylinder - (linear basis; /*/ a = 4) 



Basis Functions 

Figure 10 Lame’s cylinder - Effects of basis functions (C 1 weight 
function; lo/a = 0.6) 
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Figure 11 Rotating disk problem and the 
three models used 



Distance between nodes in the r - direction 


Figure 12 Rotating disk - Effect of nodal refinement - (C 1 

weight function; Cubic basis function; l 0 !a= 0.125) 
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Basis functions 

Figure 13 Rotating disk - Effect of basis functions - (C 1 
weight function; l 0 la = 0.125; Model-2) 



Support domain radius, {lola) 

Figure 14 Rotating disk - Effect of support domain radius - (C 1 
weight function; Cubic basis function, Model-2) 
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